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1 .  INTRODUCTION 


1.1  Purpose 

This  technical  memorandtm  discusses  mathematical  programming 
formulations  for  the  problem  of  optimizing  the  stationing  and  vectoring 
of  Combat  Air  Patrol  (CAP)  aircraft  employed  in  the  air  defense  of  naval 
task  forces. 

1 . 2  Background 

Previous  analyses  (see  references  (a)  and  (b))  are  based  on  certain 
simplifying  assumptions.  For  example,  they  assume  that  CAP  aircraft 
always  fly  lead  collision  courses  with  the  target.  There  is  no  assurance 
that  this  tactic  is  always  optimal  for  CAP  aircraft  carrying  long-range 
air-to-air  missiles  (AAMs) .  This  document  suggests  that  some  of  the 
assumptions  made  in  previous  studies  can  be  weakened  considerably  if  the 
problem  is  formulated  as  a  mathematical  program. 

1 . 3  Scope 

In  this  memorandixm  only  the  case  where  all  aircraft  and  missiles 
move  in  straight  lines  and  at  constant  speed  is  considered.  It  is 
assumed  that  the  threat  moves  directly  toward  the  center  of  the  Task 
Force  (TF) ,  and  that  CAP  aircraft  may  or  may  not  be  vectored  to  the 
target  by  friendly  airborne  or  other  long-range  surveillance  and  air 
defense  early  warning  systems.  All  action  is  treated  as  taking  place  in 
a  Euclidean  plane. 

1 . 4  Approach 

Consider  only  one  type  of  scenario,  which  is  based  on  the  disposition 
of  forces  shown  in  Figure  1.  In  this  scenario,  a  CAP  aircraft  is  stationed 
at  a  distance  from  TF  center  given  by  the  CAP  station  radius.  A  threat 
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aircraft  group  heads  directly  for  TF  center  along  a  radial  making  an 
angle  alpha  with  the  CAP  station' s  radial.  At  some  point  in  its  ap¬ 
proach,  the  threat  is  detected  and  (after  some  delay  for  threat  analysis 
and  command  decision)  the  CAP  aircraft  are  vectored  to  engage  the 
threat.  After  they  come  within  range  of  the  threat  aircraft  group,  the 
CAP  aircraft  launch  their  AAMs.  The  AAMs  intercept  the  threat  at  some 
Intercept  Range  as  illustrated  in  Figure  1.  To  be  considered  successful, 
the  intercept  must  occur  before  the  threat  aircraft  reach  their  anti¬ 
ship  missile  (ASM)  launch  range  (i.e. ,  the  Intercept  Range  must  exceed 
the  ASM  Release  Range  shown  in  Figure  1) . 

There  are  a  number  of  limitations  on  CAP  and  AAM  capability  that 
must  be  honored.  For  example,  it  is  not  practical  to  station  CAP 
aircraft  too  far  from  TF  center,  AAMs  may  have  to  be  launched  at  certain 
target  aspect  angles  (or  so  as  to  attain  certain  aspect  angles  at 
intercept) ,  and  so  forth.  Section  3  gives  a  detailed  discussion  and 
mathematical  development  of  the  constraints  considered  in  this  paper. 

The  maximum  sector  coverage  provided  by  a  CAP  stationing  and 
vectoring  policy  can  be  identified  with  the  maximum  angle  alpha  (shown 
in  Figtire  1)  that  gives  an  intercept  before  the  threat  aircraft  reach 
their  ASM  Release  Range  and  is  compatible  with  all  the  constraints  on 
CAP  and  AAM  capabilities.  Other  measures  of  effectiveness  are  con¬ 
ceivable  and  some  of  than  are  discussed  in  Section  4. 

1.5  Results 

The  results  are  a  mathematical  programming  formulation  for  a  feumily 
of  CAP  stationing  and  intercept  vectoring  problems.  A  large  variety  of 
particular  cases  can  be  obtained  by  relaxing  selected  constraints. 
Unfortunately,  the  resulting  mathanatical  programming  problems  are 
nonlinear,  and  so  can  be  solved  only  by  nvimerical  methods.  However, 
computer  routines  for  solving  such  mathematical  problems  are  available 
and  could  be  applied  to  CAP  stationing  and  vectoring  problems  of  the 
type  developed  here  (see,  for  example,  reference  (c)).  Some  suggestions 
on  ways  to  formulate  the  problems  to  facilitate  their  solution  are 
offered  in  Section  5. 
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2.  SCENARIO  AND  DEFINITION  OF  VARIABLES 


2.1  Scenario  . 

The  reader  is  referred  to  Figure  2  for  a  sketch  illustrating  the  * 

relations  among  the  key  variables  selected  for  describing  the  problem  in 
mathematical  form.  In  this  figure,  O  is  the  TF  center  and  S  is  the  CAP 
stationing  point.  The  threat  approaches  TF  center  along  the  radial  0-T 
and  at  a  speed  VT.  TD  is  the  distance  of  the  threat  fron  TF  center  when 
it  is  first  detected  by  elements  of  the  TF  -  which  may  or  may  not  be 
the  CAP  aircraft  thenselves.  RhoD  is  the  distance  of  the  threat  frati 
the  CAP  when  the  threat  is  detected,  and  omegaD  is  the  target  aspect 
angle  of  the  threat  as  viewed  from  the  CAP  at  that  time.  Lower-case  h 
is  time  that  elapses  between  detection  of  the  threat  and  commitment  of 
the  CAP  to  attack  it.  TC  is  the  distance  of  the  threat  frcsn  TF  center 
when  the  CAP  are  committed;  rhoC  and  omegaC  are  the  distance  of  the 
threat  from  the  CAP  and  the  aspect  angle  of  the  threat  as  viewed  from 
the  CAP  at  that  time.  L  is  the  point  at  which  AAMs  are  launched  from 
the  CAP  aircraft,  and  rhoF  is  the  distance  of  that  point  from  the  CAP 
station  point.  OmegaF  is  the  angle  as  shown  on  Figure  2.  VP  is  the 
speed  (assumed  constant)  of  the  CAP  aircraft  while  flying  to  the  AAM 
launch  point  and  tF  is  the  time  required  for  that  flight.  RhoL  is  the 
distance  of  the  CAP  aircraft  from  the  threat  at  the  time  the  AAMs  are 
launched  and  omegaL  is  the  aspect  angle  of  the  threat  as  viewed  from  the 
CAP  aircraft  at  AAM  launch.  TM  is  the  distance  of  the  threat  from  TF 
center  when  intercept  occurs.  OmegaM  is  the  aspect  angle  of  the  threat 
as  viewed  from  the  missile  at  that  time.  VM  is  the  speed  of  the  missile 
(assumed  constant)  during  its  flight,  and  LM  is  its  time  of  flight  from 
launch  to  intercept.  Angles  are  taken  to  be  positive  in  the  directions 
indicated  in  Figure  2.  Angles  omegaL,  omegaM  and  omegaF  may  vary  from 
-180“  to  +180°.  The  other  angles  (alpha,  omegaD  and  omegaC)  may  vary 

from  0°  to  +180°. 
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2.2  Relations  Among  the  Variables 

Many  of  the  relations  among  the  variables  presented  in  Table  1  are 
easily  derived  from  trigonometric  formulas  for  right-angled  triangles. 
The  others  follow  from  the  assumed  constant  speed  of  the  threat,  the  CAP 
aircraft,  and  the  AAMs* 

2. 2.1  Discussion  of  the  Relations  Among  the  Variables 

one  of  the  relations  C3) ,  (4)  and  (5)  shown  in  Table  1  is  redundant 
and  can  be  eliminated. 

This  memorandum  considers  only  those  cases  in  which  the  following 
parameters  have  known  fixed  non-negative  values  given  by  the  conditions 
of  the  problem;  VF,  VM,  and  VT. 

The  parameters  shown  in  Table  2  are  considered  to  be  unknowns 
controllable  by  the  TF's  anti-air  warfare  system,  svibject  to  the  re¬ 
lations  among  them  given  in  Table  1  and  to  appropriate  constraints  to  be 
discussed  in  Section  3.  Because  of  the  relations  among  them  and  to  the 
constraints,  not  all  of  the  controllable  variables  can  be  adjusted 
independently . 
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TABLE  1.  RELATIONS  AMONG  THE  VARIABLES 


Pg  sin  ® 


Pq  sinw^ 


P_  sin  ® 
S 


p  ^  sinuj^  =  P  jj  sino) 


p  sin  0)  ■+  P  T  “  T  “  ^  s  ^ 


Pp  sin  Up  +pj^  sinUj^  =pg  sin  a 


Pd  cos  a)p  =  Tp  -  Pg  cos  a 


p^  cos  =  T^  -  p  g  cos  a 


p^  cos  ^  ~  "^L  ”  ^S 

Pm  “m  Pp  “f  “  '^M  ■  ^s  “ 


pp  cos  Wp  =  Tp  -  Pg  cos  a 


T„  =  -  v„  h 

C  D  T 


T  =  T  -  V  t 
L  C  T  F 


T  =  T  -  V  t 
L  T  M 


’’f  ■  •  I’m  “m 


p  =  V  t 
F  F 


p  =  V  t 
'H  H 
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3.  CONSTRAINTS 


3.1  Forced  Constraints 

The  forced  constraints  that  seem  most  reasonable  to  impose  are 
shown  in  TeOble  3.  Others  may  be  appropriate  in  different  contexts.  Not 
all  of  them  will  actually  be  effective  in  any  given  problem  situation. 
Those  that  are  not  to  be  effective  can  either  be  omitted,  or  else  their 
bounds  can  be  taken  so  large  Cor  so  small)  that  the  constraint  has  no 
effect  on  the  solution. 


One  interpretation  of  the  forced  constraints  will  be  offered  here. 
Others  are  possible  and  may  be  useful,  depending  on  circumstances.  RD 
can  be  viewed  as  the  maximum  detection  range  of  the  threat  by  CAP  air¬ 
craft,  and  RC  as  a  maximum  range  for  identifying  the  threat  and  estimating 
its  raid  size  by  the  CAP  aircraft.  RL  can  be  considered  a  maximum  AAM 
launch  range,  RM  a  maximum  AAM  flight  range,  RF  a  maximum  CAP  flight 
distance,  and  RS  a  maximum  CAP  station  radius. 


YD,  YC,  YL  and  YM  are  the  maximum  distances  of  the  threat  from  TF 
center  at  which  it  can  be  detected,  CAP  can  be  committed  to  engage  it, 
AAMs  can  be  launched,  and  AAM  intercepts  can  occirr  (respectively).  YF 
is  a  constraint  that  has  been  inserted  as  analogous  to  those  just 
mentioned,  although  no  physical  significance  has  yet  been  attached  to 
it.  YASM  is  the  threat's  ASM  Release  Range.  H  is  the  minimum  time  for 
analyzing  the  threat  and  deciding  to  commit  CAP  to  engage  it.  Con¬ 
straints  (14)  through  (17)  are  restrictions  on  the  target  aspect  angles 
at  detection,  CAP  commitment,  AAM  launch,  and  intercept  (respectively) . 
In  each  case,  they  prevent  these  angles  from  being  too  large.  For 
example,  constraint  (14)  requires  that 


-  0)  < 
DO  - 


<  CO  # 
-  DO' 


and  so  forth. 
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TABLE  3.  FORCED  CONSTRAINTS 


(1) 

"c  -"c 

C2) 

(3) 

C4) 

Pp  — 

(5) 

Ps  -“s 

(6) 

D  —  D 

C71 

^c-^^c 

C8) 

C9) 

T  <  Y 

M  -  M 

CIO) 

T  >  Y 

M  -  ASM 

Cll) 

T„  <  Y^ 

F  —  F 

(.12) 

h  ^  H 

(13) 

cos  ^  cos 

(14) 

cos  (0^  ^  cos 

(15) 

cos  0)^  ^  cos  (O^Q 

(16) 

cos  0)^  >  cos  a)j^Q 

(17) 

cos  COp  ^  cos  UpQ 

(18) 
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Constraints  (14).  through  U?)  are  couched  in  terms  of  limitations 
on  the  cosines  of  the  angles,  rather  than  in  terms  of  the  angles  them¬ 
selves.  The  reason  for  this  will  become  more  apparent  in  Section  5, 
where  it  is  suggested  that  algebraic  functions  be  substituted  for 
trisonometric  ones  via  the  substitutions : 


cos  (angle)  =  x 
sin  (angle)  =  y 


and  the  introduction  of  the  following  relation  among  these  new  variables 


2 

+  y 


=  1 


3.2  Natxiral  Constraints 

Those  natural  constraints  that  seam  most  reasonable  to  impose  are 
shown  in  Table  4.  Natural  constraints  (7)  through  (12)  are  stated 
explicitly  in  anticipation  of  the  suggested  change  from  trigonometric 
to  algebraic  functions  mentioned  above  and  discussed  in  Section  5. 

3.3  Discussion  of  Constraints 


Natural  constraints  (4)  and  (4a)  are  equivalent  because  VM  is  non. 
negative.  Similarly,  natural  constraints  (5)  and  (5a)  are  equivalent 
because  VF  is  non- negative. 


Forced  constraints  may  be  dropped  qualitatively  by  omitting  a 
constraint,  and  this  is  sometimes  the  desirable  way  to  proceed.  How¬ 
ever,  it  is  also  sometimes  desirable  to  drop  a  constraint  quantitatively 
by  setting  its  bound  to  such  a  large  (or  small)  value  that  the  constraint 
has  no  effect  on  the  solution. 

This  memorandum  considers  only  those  cases  in  which  the  bounding 
values  shown  in  Table  5  have  fixed  values  given  by  the  conditions  of 
the  problem.  Except  for  the  trigonometric  functions,  these  values  will 
be  non-negative  for  most  problems. 
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TABLE  4.  NATURAL  CONSTRAINTS 


Pd  -  ° 

Pc  1  0 

Pl  >  0 

Pm  1  0  1  0^ 

Pp  >  0  1  0)- 

Pg  >  0 


2  .2 

cos  0)^  +  sin  =  1 


2  2 
cos  0)^  +  sin  (A)^  =  1 

c  c 


2  2 
cos  0^  +  sin  0^  =  1 


2  2 
cos  (o„  +  sin  a)„  =  1 
M  M 


2^-2  , 
cos  a)„  +  sin  a)„  -  1 
F  F 


2^.2  . 
cos  a  +  sin  a  =  1 


0®  <  0)  <  180° 

—  D  — 

0°  <  w  <  180° 

—  C  — 

-180°  <ca^  <180° 
—  L  ““ 

-180“  <  0).,  <  180“ 
—  M  — 

-180“  <  a)„  <  180“ 
—  F  — 

0“  <  a  <  180“ 


(1) 

t2) 

(3) 

C4) C4a) 

(5)  (5a) 

(6) 

(7) 

t8) 

C9) 

(10) 

(11) 

(12) 


» 
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TABLE  5.  LIMITING  VALUES 


H 

COS  0) 

cos  0) 

( 

cos  0) 
cos  0)^ 
cos  0). 


DO 


CO 


LO 


MO 


FO 


ASM 


4 


4 
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4.  OBJECTIVE  FUNCTIONS 

Table  6  shows  some  possible  objectives  to  be  achieved  by  manipu¬ 
lating  the  controllable  variables.  A  choice  of  one  of  them  (or  of 
anotiier  not  shown  in  Table  6)  depends  on  a  mixture  of  considerations 
involving  mathamatical  convenience  and  appropriateness  to  the  applica¬ 
tions  being  contemplated  for  the  solution. 

Objectives  (1)  and  (21  are  equivalent,  and  both  ask  for  the  largest 
angle  subtended  by  the  region  defended  by  the  CAP  station.  Objective  (31 
asks  for  the  maximxan  arc  length,  and  objective  (41  for  the  maximum  length 
of  "front,”  defended  by  the  CAP  station.  Objectives  (51  and  (61  ask  for 
the  maximum-area  circular  sector  and  triangle  defended  by  the  CAP  station. 
Objective  (71  asks  for  the  maximum  time  delay  attainable  —  if  this  objec¬ 
tive  is  adopted  it  is  advisable  to  add  a  forced  constraint  preventing 
alpha  from  being  too  small,  since  otherwise  the  maximum  h  will  probably 
occur  at  alpha  =  0.  Such  a  bound  on  alpha  is  also  advisable  if  objective 
(81,  i.e.,  maximizing  the  distance  between  the  intercept  point  and  TF 
center,  is  selected.  Objectives  (91,  (101  and  (111  are  variations  of 
objective  (8}  that  autcmatically  prevent  the  optimizing  alpha  from 
becoming  zero.  If  objective  (121  is  adopted,  then  it  may  be  advisable 
to  add  a  forced  constraint  that  prevents  rhoS  from  becoming  too  small, 
since  otherwise  the  minimum  may  occur  at  rhoS  =  0. 
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TABLE  6.  CANDIDATE  OBJECTIVES 


Maximize  a 
Minimize  cos  a 
Maximize  ap 

o 

Maximize  sin  a 

Maximize  ap^ 

2  12 

Maximize  p  sin  d  cos  =  Pe  sin  C2a) 

Maximize  h 

Maximize 

M 

Maximize  Y„  a 
M 

.  .  2 

Maximize  Y„  a 
M 

.  .  2 

Maximize  Y^  a 

Minimize  cos  a 


ti) 

C2) 

t3) 

C4) 

C5) 

(6) 

(7) 

C8) 

19) 

(10) 

(11) 

(12) 


* 


* 
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5.  DISCUSSION 

If  the  CAP  aircraft  are  not  equipped  with  long-range  AAMs,  then  ^ 

take  RM  =  0.  ^  forced  constraint  C4) ,  this  will  immediately  imply 

rhoM  =  0  and  hence  by  relation  a6).  of  Table  1  tM  =  0.  Then  it  follows  . 

from  relations  a3)  through  aS).  that  TL  =  TM  =  TF .  And  then 
rhoL*sin  ComegaL).  =  rhoL*c6s  Comeg aL).  =  0;  hence  rhoL  —  0. 

It  is  strongly  suggested  that  before  a  nonlinear  optimization  code 
is  applied  to  this  problem,  all  trigonometric  functions  be  replaced  by 
their  component  coordinate  values  and  appropriate  additional  constraints 
be  introduced  to  express  the  relation  among  the  component  coordinate 
values.  Thus,  for  example,  we  \arge  the  replacements; 

cos  (alpha)  =  AX, 
sin  (alpha).  =  AY, 

AX^  +  AY^  =  1, 

and  similarly  for  all  the  other  angles.  The  aim  here  is  to  rid  the 
problem  of  trigonometric  functions,  replacing  them  with  expressions  that 
are  no  worse  than  quadratic  in  the  variables.  When  this  is  done,  all 
the  necessary  numerical  evaluations  of  functions  and  their  gradients  can 
be  done  by  multiplication  and  addition,  which  are  performed  much  quicker 
than  the  evaluation  of  trigonometric  functions. 


It  is  sometimes  hard  to  find  an  initial  feasible  solution,  espe¬ 
cially  to  nonlinear  problems  such  as  this  one.  Nevertheless,  in  this 
case,  an  initial  feasible  solution  can  conveniently  be  found  by  the 
following  argument.  First,  it  seems  intuitively  obvious  that  there  is  a 
feasible  solution  for  the  original  nonlinear  problem  if  and  only  if 
there  is  a  feasible  solution  for  the  "collinearized"  problem  obtained  by 
setting  all  of  the  angle  variables  to  zero.  It  is  easily  seen  that  all 
the  relations  among  variables  and  constraints  for  such  a  "collinearized 
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problem  are  either  linear  equations  or  linear  inequalities.  So  if  we 
adopt  them  as  the  constraints  and  ask  for  a  maximum  value  of  rhoS,  we 
find  (see  Table  7)  that  we  have  a  linear  programming  problem,  the 
solution  of  which  (if  there  is  one)  provides  an  initial  feasible  solu¬ 
tion  for  the  original  nonlinear  programming  problem. 

Various  search,  quasilinearization,  or  feasible  direction  methods 
can  be  used  to  improve  on  any  current  feasible  (but  nonoptimal)  solution. 
See  reference  (d)  for  a  survey  of  potentially  applicable  techniques  and 
references  to  other  sources.  If  this  work  is  to  be  pursued,  it  is 
recommended  that  either  a  Lagrange  multiplier  or  feasible  direction 
method  be  tried. 


* 
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TABLE  7.  COLLINEARIZED  PROGRAMMING  PROBLEM 


Maximize:  P, 


Subject  to; 


p  =  T  - 
D  D  S 

Pc  =  Tc  -  Pg 

'^L  *  -  ■'i,  -  "s 

‘'m  ^  '  ’’m  -  "s 

-  Pc 

F  F  S 

=  Tc  -  v„  h 

C  D  T 

T  =  T  -  V  t 
L  C  T  F 

T  =  T  -  V  t 
M  L  T  M 

"^F  “  "^M  ”  ^M 

P-  =  v_  tp 


F  F 
5 

M 


•  ’h  'm 


Table  3,  Forced  Constraints  (1)  through  (.13)  apply. 
Table  4,  Natural  Constraints  (1)  through  C6)  apply. 
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